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Abstract 



Effective field theory, including pions, provides a consistent and systematic 
description of nucleon-nucleon strong interactions up to center-of-mass mo- 
mentum p ~ 300 MeV per nucleon. We describe the inclusion of hadronic 
parity violation into this effective field theory and find an analytic form for 
the deuteron anapole moment at leading order in the expansion. 



I. INTRODUCTION 

An increasing amount of both experimental and theoretical interest is being focused on 
nuclear parity violation. The problems in reproducing parity violating observables in light 
nuclei (such as ^'^F [0]) and in the nuclear anapole measurement of cesium |^ by parameters in 
single boson exchange potentials suggest that we re-examine the theoretical framework 
with which hadronic parity violation is studied. The low momentum transfers involved in 
these parity violating processes makes them amenable to an effective field theory treatment. 
After much effort |5|^3|, a consistent power counting and calculational scheme has been 
developed which provides a systematic treatment of the low-momentum nucleon-nucleon 
strong interactions p3[, up to a center-of-mass momentum for each nucleon of p ~ 300 MeV. 



Analytic expressions for A^A^ scattering phase shifts and deuteron observables P^,p5[ have 



been obtained which describe the data well. Using this technique we construct an effective 
field theory describing weak interactions in the two-nucleon sector, and as an example obtain 
an analytic expression for the deuteron anapole moment at leading order in the expansion. 



*This paper is published in Nucl. Phys. A644 235-244 (1998). The last section contains an 
erratum and addendum to the published version. 

^On leave from the Department of Physics, Duke University, Durham NC 27708. 
rpsSredhook . phys . Washington . edu . 
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II. STRONG INTERACTION AMPLITUDES 



The strong interactions of the pions are described by the Lagrange density 



£2 £2 
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where the pion fields are incorporated in a special unitary matrix, 
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with / = 132 MeV. The strong interactions of the nucleons are described by a Lagrange 
density that we write as 



where contains n-body nucleon operators. 

The one-body terms in the Lagrange density are 
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where qa = +1.25 and the ellipses denote operators with more insertions of the light quark 
mass matrix and more spatial derivatives. The two-body Lagrange density contributing to 
S-wave interactions may be written as 
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where the ellipses denote operators involving more insertions of the light quark mass matrix, 
meson fields, and spatial derivatives. The Pj^o and Po^a are spin-isospin projectors defined 
by 



Pi,0 = 



Po,a = —?=(y2T2Ta 



1 

— ( 
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(2.6) 



The coefficients have been determined from fits to the S-wave phase shifts in both the ^5*0 
and ^5*1 channels ||2^ . Renormalized at /i = in the power divergence subtraction scheme 
(PDS) it is found that 
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Cj'^°) = -3.34 fm^ , = -0.42 fm^ , cf^°^ = 3.24 fm^ , 

C^'^i) = -5.51 W , Df^'^ = 1.32 W , Cf'^'^ = 9.91 fm^ . (2.7) 

The relative importance of an operator to a particular observable is given by the running 
of the operator in the infrared. Dimensional regularization with the PDS scheme provides 
a simple and consistent determination of this running. Consider a general operator arising 
from strong interactions that connects states with angular momentum L and L' . If this 
operator has m insertions of the light quark mass matrix and 2d = L + L' + 2n spatial 
gradients, we will denote its coefficient by C^;f ^. It was shown in p3| , p4| that such a 



coefficient scales like /x for the ^5*0 and ^5*1 — ^Di partial waves, and like ii^ for 

all other partial waves. For nucleon nucleon scattering in the ^5*0 channel the momentum 
independent contact term with no insertions of niq has coefficient Cq'^"'' = C'^q'' ~ and 
is the leading contribution. Each loop graph in the time-ordered product of two insertions 
of this operator scales like ixC'^q^ ~ 1. Therefore the bubble chain of such insertions 
constitutes the leading amplitude. At subleading order two additional operators appear. 
There is a contribution from a single insertion of an operator with (i = 1, m = (coefficient 
C2 ° ), and another from a single insertion of an operator with ci = 0, m = 1 (coefficient 

°'). At the same order, there are contributions from the exchange of a single potential 
pion, which scales like For the remaining discussion we will represent the small expansion 
parameters, the external momentum p, the quark masses mg, and the renormalization scale 
/i, by the generic quantity Q. 



III. PARITY VIOLATING INTERACTIONS 

The standard model of electroweak interactions gives rise to parity violating four-quark 
operators that exhibit enhanced symmetry in particular limits. For instance, in the limit that 
the angles of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and the weak mixing angle 
vanish, and the mass of the two quarks in each generation are degenerate, the standard 
model Lagrange density has a global SU{2)l ® SU{2)fi chiral symmetry [^. All parity 



violation in this limit must be A J = 0. The large difference between the mass of the charm 
and strange quarks breaks this symmetry. For hadronic processes the strange quarks can be 
considered dynamical while the charm quarks are "integrated out" of the theory. Therefore, 
in the limit of vanishing mixing angles in the CKM matrix and in the neutral gauge boson 
sector, all AI = 1 parity violation comes from strange quarks. In the real world with non- 
zero mixing angles, the up and down quarks also contribute, but with coefficients suppressed 
by sin^ 9^, and by the Cabibbo angle. 

Excellent reviews of parity violating observables in the context of meson exchange models 
of weak nucleon-nucleon scattering appear in p6|-pO[] . In this work we do not reproduce the 



results obtained in potential models for various parity violating observables, but instead 
focus on the effective field theory description of parity violation. This will facilitate the 
systematic inclusion of higher order effects such as relativistic corrections and dynamical 



meson exchange. A first step was undertaken in |^ where a Lagrange density consistent 
with chiral symmetry was constructed for the single nucleon sector. As with any effective 
field theory constructed to reproduce the low-energy behavior of some unknown or unsolved 
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theory, there are couphngs in the theory that must be determined by experiment in order 
for the theory to be predictive. 

The parity violating four-quark operators generated by the standard model of electroweak 
interactions can be classified according to how they transform under SU {2) SU {2) chiral 
transformations By explicit construction they have decompositions (1, 1), (3, 1) © (1, 3), 
and (5, 1) © (1, 5), and so can be unambiguosly classified by their transformation under the 
isospin subgroup. Operators transforming as (3, 3) also occur but do not violate parity. The 
following operators |31[ 
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(3.1) 



transform as X — > UXW under chiral transformations. Rewriting the Lagrange density of 
| 3l| in terms of two-component non-relativistic nucleon fields, so that they are expressed in 
terms of the same field operators as the strong NN Lagrange density in eqs.( |2.4|) and ( p. 51) , 
we haveQ 
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where the ellipses represent higher order terms in the chiral and momentum expansion. A^ 
is the axial vector meson field A^ = ^ {^^D^^ — ^-D/^^^), with = df^ + ieA^, and Tab is 
defined in 113111 to be 
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1 

Vo 








(3.3) 



In the power counting we employ, the size of the coefficients in the Lagrange density of 
eq. ( W^ i^wNN^ hy'^''^\ and /i^'^"*) is set by naive dimensional analysis (NDA). Were the co- 
efficients to dramatically deviate from their NDA size the power counting must be modified. 
Expanding the interactions in eq. (|3.2|) to 0{n) gives 



'-ttNN 



-i h^^lj^ 7r+ p'^n + h.c 



(3.4) 



where the hy-type interactions involving a single vr do not contribute due to the conservation 
of the vector current. The strong interaction between two nucleons induced by the exchange 



^ Note that the sign of h^^J^j^ is opposite that used in |31] 



4 



of a single potential pion scales like Q*^ in the effective field theory; two factors of momentum 
from the derivative coupling and two powers of q and in the pion propagator give 

X (ff ■q)(ff -q) , . 

V{q) ~ 5 5 — . (3.5) 

The weak interaction induced between two nucleons from the single insertion of the oper- 
ator with coefficient h^^j^j^ scales like Q~^, since the weak ttNN vertex is independent of 
momentum, 

y(q)~^^ . (3.6) 



The meson-nucleon vertices induced by the other operators in eq.( |3.2| ) involve more deriva- 
tives so their contribution to the weak nucleon nucleon interaction is suppressed by additional 
powers of the expansion parameter Q. If for some reason /iItatat is much smaller than indi- 
cated by NDA, as suggested by the ^^F experiments while the others coefficients remain 
their NDA size, then this weak one pion exchange interaction will no longer be the leading 
order. 

In addition to the contribution from the exchange of one or more mesons, there are con- 
tributions to the weak interaction between nucleons coming from weak four-nucleon contact 
operators. The leading order weak interaction (CP invariant) between nucleons in an S-wave 
in both initial and final states is described by the Lagrange density 

^ = vis',^So;0^ {N^Po,aNy N^{Po,a,Xl-Xl}TN + h.C. 

^ Visl^So-fl {p^ (^2n + a2pf {p^ 02P T^" - (T2n -K^^ + h.c. + 0{tx^) . (3.7) 



Notice that this does not generate an interaction involving only two nucleons, but does 
give rise to an interaction between two nucleons and pions. We have defined the modified 
anti-commutator by 

{A, B}t = AB + A . (3.8) 

Coefficients of these operators are defined to have the form (2s, +11 , where 

the nucleon quantum numbers are defined to be those in the absence of pions, and n is 
the number of and rrig insertions. The /3-function for these operators shows that that 
coefficient rji^^^-^g^.Q scales like In general, operators in this channel with n insertions of 
and m insertions of the light quark mass matrix scales like ^-2('T-+ni+i)_ p^-^ next order 
there are contributions from operators involving more insertions of the light quark mass 
matrix, insertions of the axial vector meson field and more derivatives acting on the 



nucleon fields. This is similar to those found in eq.( |2.5| ). There is no analogous interaction 
involving nucleons in the ^5*1 channel that is CP invariant. 

The dominant contribution to low-energy parity violation in pp scattering is from the 
^5*0 — ^Pq weak amplitude interfering with the strong scattering amplitude. The four nucleon 
operators for such weak scattering are, for AI = 



5 



5 P, , 



+ h.c. 



and for AI = 1 



(3.9) 



P,; - Id P,; 



N + h.c. 



N + h.c. . (3.10) 



We have not shown the AI = 2 Lagrange density. The coefficients of these operators, Visl^Po O^ 
Vist^Pg-Qj s-^d rjig^sp^.Q scale hke and hence contribute to the pp scattering amphtude at 
order Q^. In general, operators in this channel with n insertions of and m insertions of the 
light quark mass matrix scales like At next order in the expansion terms involving 

more insertions of spatial gradients and the light quark mass matrix give a contribution at 
order Q^, suppressed by only one power of Q in the expansion. Therefore, at next-to- 
leading order (NLO), there are contributions from both subleading terms in the strong 
interaction and from subleading terms in the weak interaction. The construction of the 
four nucleon operators contributing to weak processes between higher partial waves follows 
straightforwardly from the above discussions, and we will not detail them in this work. 



IV. THE ANAPOLE MOMENT OF THE DEUTERON 

As an example of a parity violating quantity that can be described using this effective 
theory, we will calculate the anapole moment of the deuteron at leading order in the expan- 
sion. This is of interest because of the possibility of measuring or constraining the coefficient 
/i^^jvAf- Attempts to determine h^^jfj^ through parity violating observables in nuclei have been 
frustrated by the uncertainties inherent in complicated nuclear systems. Parity violating ob- 
servables in the interaction between leptons and single nucleons also receive contributions 
from /ilrATAT but there are contributions from both isoscalar and isovector mediated in- 
teractions as well. In contrast, the deuteron has only isoscalar interactions, dominated by 
the strangeness content of the deuteron. Fortunately, there are currently two experiments 
under consideration p3| , p4| which will look at the theoretically clean deuteron system. The 



rip — > d'-y process 33 1 has been treated in this effective field theory in ref. [pq] and parity 



violation in electron-deuteron scattering is the subject of this section. 

The anapole moment of a particle is its spin dependent coupling to the electromagnetic 
field. For the nucleon 

C = a'' N d^F'^'^ . (4.1) 

The equations of motion for the electromagnetic field, d^F^^ = e\E'7'^\E', allow this operator 
to be written entirely in terms of local four- Fermi contact terms. It is convenient to keep the 
operator in the form of eq. (|4.1|) for calculational purposes. Parity violating lepton hadron 
scattering that depends upon the spin of the hadronic target receives contributions from 
both the anapole moment of the object and from the exchange of gauge bosons between 
the lepton current and the hadronic target. The anapole moment of elementary particles 
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FIG. 1. Leading order diagrams contributing to the deuteron anapole moment. The crossed 
circles denote operators that create or annihilate two nucleons with the quantum numbers of the 
deuteron. The solid square denotes the weak operator with coefficient h^j^j^ and the solid circle 
denotes minimal coupling to the electromagnetic field. Wavy lines are photons, solid lines are 
nucleons, and dashed lines are mesons. 



depends upon the choice of gauge |36|, since its contribution cannot be separated from 
the contributions of IV-box and Z-loop graphs, which all occur at order g"^ in the electroweak 
coupling. However, since the anapole moment of the hadrons (composite objects) is dom- 
inated by the long distance behavior of QCD, and the weak vertices at the hadronic level 
are independent of the choice of gauge, these are well-defined contributions to the spin 
dependent weak scattering amplitude. 

The anapole moment of the deuteron is Ad, the coefficient in 

^ = '"^0^ eabcD''^D''d,F^' , (4.2) 

where ieabcD°''^D^ is the deuteron spin operator. The leading contribution to Ar, comes from 
weak meson exchange starting at order Q~^, dominated by the operator with coefficient 
( see fig. (|I|).) The weak four nucleon operators start contributing at order and we 
will neglect them in this work. 

Explicit computation of the graphs shown in fig. (|l]) leads to a deuteron anapole moment 

of 



- TTT-f 



6(m, + 27)2 



(4.3) 



where 7 = y/M^B is the binding momentum of the deuteron with B the deuteron binding 
energy. The first term in eq. ( |4.3| ) is from the graphs that give the dominant contribution 
to the nucleon anapole moment (the last four graphs in fig. (|l])), and the second term 
is the contribution from potential pion exchanges between the two nucleons in the deuteron 
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(the first six graphs in fig. (||)). These two terms are the same order in the expansion, 
but the contribution from the individual nucleon anapole moments is numerically larger. 
Higher order corrections to this result that are parametrically suppressed by the expansion 
parameter(s) arise from higher dimension operators in both the weak sector and in the strong 
sector. These corrections are typically on the order of 30%. In addition to the contributions 
from the higher dimension weak meson operators in ( p.2|) and the contributions from the 
weak four nucleon operators there is a contribution from four nucleon anapole operators, 
e.g. 

£ = ie,,u{N^P,,^N)\N^P,^,N) d^F^' + ... , (4.4) 

that contributes at next-to-next-to-leading order in the Q expansion. The coefficient A^j^ 
scales like in the PDS scheme. 

The deuteron spin dependent amplitude for electron deuteron scattering receives contri- 
butions from both the deuteron anapole moment and from the direct exchange of a Z° (and 
the associated radiative corrections). The spin-dependent matrix element for eD scattering, 
at tree-level in the standard model is 

A^'""^ = 10-^ [^Ag^s - 0.17 (^^j j ^el^e e^'Ke,. , (4.5) 

where e* and are the polarization vectors of the final and initial deuterons and g^s is 
the strange axial matrix element of the nucleon. Loop corrections in the standard model 
are fractionally larger than one would naively expect largely because of the suppression of 
the tree-level amplitude by a factor of (1 — 4sin^6'u,) The estimated value of /ilrjvAr is 



in the range — > 11 x 10 ^ from the work of DDH [^, while naive dimensional analysis 
estimates give h^^lfj^ ~ 5 x 10^'' ||31|j3^. Present determinations of g^s from lepton scattering 



38| , |39[| , including SU(3) breaking effects, are ~ —0.15 ± 0.15 (e.g. [0). This suggests that 
the deuteron anapole moment contribution to A^^^^ is comparable to the contribution from 
Z° exchange. However, for values of h^^li^ suggested in or larger, the deuteron anapole 
moment makes the larger contribution to A^'"\ While it is clear that this will not provide 
the cleanest extraction of h^^^li^, we encourage the experimental community to continue 
their investigation of the parity violating spin dependent interactions of the deuteron (e.g., 
p^ ) in order to constrain h^^j^j^ and help resolve the current controversy surrounding this 



parameter. (Of course, if h^^^li^ turns out to be much smaller than NDA suggests, then the 
hadronic contribution we have computed here will not be the dominant one.) 



V. CONCLUSIONS 

We have shown how to systematically include hadronic parity violation from the standard 
model of electroweak interactions into an effective field theory description of nucleon nucleon 
interactions, including photons and pions. Exchange of a single potential pion with coupling 
h^T^lij^ is the leading contribution to the weak scattering of two nucleons in the effective field 
theory expansion. The most general parity violating four nucleon operators categorized by 
partial wave and isospin change were discussed and the lowest dimension operators presented. 
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The leading four nucleon operator that violates parity has at least one pion in addition to 
nucleons. As an example we computed the leading contributions to the deuteron anapole 
moment, which could be used to place a constraint on /i|rjvAr- 



We would like to thank Ernest Henley and Michael Ramsey-Musolf for useful discussions. 
This work is supported in part by the U.S. Dept. of Energy under Grants No. DE-FG03- 
97ER4014 and DE-FG02-96ER40945. 
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VI. ERRATUM AND ADDENDUM 



Recent work by Khriplovich and Korkin (KK) |^ has demonstrated that a leading 
order (LO) contribution to the deuteron anapole moment was omitted in our paper . 
In addition to this omission there is a factor of —2 error in one of the contributions. In 
this "Erratum and Addendum" we correct the errors, and discuss the chiral limit of this 
observable in more detail. 

There is a contribution to the deuteron anapole moment arising from the isovector mag- 
netic moment of the nucleon which was omitted in 
is described by the lagrange density 



, as detailed in ^T[]. This interaction 



A, 



2M, 



N 



-N'' (ko + K1T3) a ■ BN 



(6.1) 



where kq = |(Kp + k„) and ki = ^{np — are the isoscalar and isovector nucleon magnetic 
moments in nuclear magnetons, with 



Kr, = 2.79285 



-1.91304 



(6.2) 



The magnetic field is conventionally defined B = V x A. The experimentally determined 
values of kq and ki include the pion loop corrections that we need to this order. 



Using the same definition of the deuteron anapole moment as given in |42] of 



/: = iA 



D 



our corrected result becomes 
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(6.4) 



where 7 = with B the deuteron binding energy. The first term in eq. (|6.4| ) is gener- 

ated by the nucleon isovector magnetic moment. The second term is the leading contribution 
from the nucleon anapole moment, which is purely isoscalar. The third term in eq. ( |6.4| ) 
arises from the pion exchange interaction between the nucleons. Combining the last two 
terms gives the somewhat more compact expression, 
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1277/ 



+ 7 



+ 



2m^ + 97 



(m^ + 27)2 6(m^ + 27)^ 



(6.5) 



Inserting the appropriate values for constants into eq. ( |6.5|) , the tree-level parity violat- 
ing matrix element for electron-deuteron scattering that depends upon the deuteron spin 
becomes 



ttNN_ 
-7 
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1 



Ml 



(6.6) 



This expression replaces eq. (4.5) in The coefficient of h^^l;j^ has been changed from 



0.17 to 0.21. The individual contributions to the numerical value of 0.21 that appears 
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in eq. ( p.6|) are 0.17 from the nucleon isovector magnetic moment, 0.07 from the anapole 
moment of the nucleon, and —0.03 from minimal coupling interactions between the nucleons. 
The amplitude is dominated by the interaction with the nucleon magnetic moment. It is 
important to note that the expression given in eq. (|6.6|) is valid only for momentum transfer 
|k| ^ TjW^tt- For momentum transfers larger than this, there will be an additional form 
factor that multiples the coefficient of h^^^ 



NN- 



In the limit of the deuteron having zero binding energy, 7 — >■ 0, becomes 



(1) 
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egAK 



m2 



Ki + 



(6.7) 



which agrees with the analogous limit of the expression given in |^ . The expressions given 
in 1^ are only valid in the limit that the deuteron 7 -C mir, as is stated below eq. (12) in 



41| . The approximation they make is justified for the actual numerical values of 7 and m-,r. 

A nontrivial check of our calculation can be made by considering the chiral limit <^ 
7 <^ A^. Since the deuteron is an isoscalar object, with a vanishing one-pion strong coupling, 
there cannot be a 1/m^ divergence in this limit. Here the charged pions cannot resolve the 
deuteron as a bound state of two nucleons and a logarithmic dependence on rrin is the most 
divergent behavior that can appear. For m,r ^ 7 eq. (|6.5|) becomes 



A 



D 



/\/f2 



487r/7 



3 

+ - 



(6.8) 



which has the expected form. The behavior shown in eq. ( |6.8| ) does not agree with the 
results obtained in who find 



A 



D 



1 



— \ - 

24/ 7 m 



(6.9) 



In fact, since it was assumed that 7 -C in deriving the results of [0, the -C 7 limit 
cannot be taken, and the dependence on ^jm-,^ is correct only in the limit of small 7/m,r- 

Our calculation does not include contributions beyond leading order in l/M^r or 1/A^ 
in the expressions for either the nucleon or deuteron anapole moment. Beyond the terms 



proportional to /i^^j^^ that contribute at higher orders, as discussed in there are also 



terms generated from additional weak interaction vertices involving pions and nucleons. 
A detailed discussion of such interactions and their contribution to the nucleon anapole 
moment is given in |^|. Therefore, the higher order correction included in [JT| is only one 



of several contributions, and only provides an indication of the size of higher order terms. 
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